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Abstract 

In this study, a numerical modeling approach was used to investigate the current-voltage behavior of conductive 
nanoplatelet-based nanocomposites. A three-dimensional continuum Monte Carlo model was employed to 
randomly disperse the nanoplatelets in a cubic representative volume element. A nonlinear finite element-based 
model was developed to evaluate the electrical behavior of the nanocomposite for different levels of the applied 
electric field. Also, the effect of filler loading on nonlinear conductivity behavior of nanocomposites was 
investigated. The validity of the developed model was verified through qualitative comparison of the simulation 
results with results obtained from experimental works. 
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Background 

In recent years, the nonlinear electrical conductivity be- 
havior of nanoparticle-modified polymers has received 
considerable attention by researchers, and several studies 
have been carried out to investigate the current-voltage 
characteristics of conductive nanocomposites. Even though 
several studies investigated the nonohmic conductivity 
behavior of insulator polymers filled with conductive 
spherical and stick-lilce inclusions [1-5], to the best of the 
authors' knowledge, all of the research in this field has 
been limited to experimental works. Experimental research 
devoted to the electric properties of nanoplatelet-based 
nanocomposites investigated the electrical conductivity 
of polymers with exfoliated graphite sheets with sizes 
varying from a few microns to several hundreds of a mi- 
cron [6-10], which allows for only limited prediction of 
the conductivity behavior of nanocomposites with submi- 
cron size inclusions. 

These limitations motivated the present authors to con- 
duct a numerical study to investigate the current-voltage 
behavior of polymers made electrically conductive through 
the uniform dispersion of conductive nanoplatelets. Spe- 
cifically, the nonlinear electrical characteristics of conduct- 
ive nanoplatelet-based nanocomposites were investigated 
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in the present study. Three-dimensional continuum Monte 
Carlo modeling was employed to simulate electrically con- 
ductive nanocomposites. To evaluate the electrical proper- 
ties, the conductive nanoplatelets were assumed to create 
resistor networks inside a representative volume element 
(RVE), which was modeled using a three-dimensional non- 
linear finite element approach. In this manner, the effect of 
the voltage level on the nanocomposite electrical behavior 
such as electrical resistivity was investigated. 



Methods 

Monte Carlo modeling 

Theoretically, a nanocomposite is rendered electrically 
conductive by inclusions dispersed inside the polymer 
that form a conductive path through which an electrical 
current can pass. Such a path is usually termed a percola- 
tion network. Figure 1 illustrates the conductivity mechan- 
ism of an insulator polymer made conductive through the 
formation of a percolation network. In this figure, ele- 
ments in black, white, and gray color indicate nanoplate- 
lets that are individually dispersed, belong to an electrically 
connected cluster, or form a percolation network inside 
the RVE, respectively. Quantum tunneling of electrons 
through the insulator matrix is the dominant mechanism 
in the electric behavior of conductive nanocomposites. 
Figure 2 illustrates the concept of a tunneling resistor for 
simulating electron tunneling through an insulator matrix 
and its role in the formation of a percolation network. 
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Figure 1 Schematic of a representative volume element 
illustrating nanoplatelets (black), clusters (white), and 
percolation network (gray). 



Electron tunneling through a potential barrier exhibits 
different behaviors for different voltage levels, and thus, 
the percolation behavior of a polymer reinforced by con- 
ductive particles is governed by the level of the applied 
voltage. In a low voltage range {eV~ 0), the tunneling re- 
sistivity is approximately proportional to the insulator 
thickness, that is, the tunneling resistivity shows ohmic 
behavior [11]. For higher voltages, however, the tunnel- 
ing resistance is no longer constant for a given insulator 
thickness, and it has been shown to depend on the ap- 
plied voltage level. It was derived by Simmons [11] that 



the electrical current density passing through an insula- 
tor is given by 

/ = /o |l exp (-^^') -{l + eV) exp \^-A (^-AA^ ] | 

(1) 

where 

Jo = el2nh{^Asf and A = {4n^As/h){2m)^ 

Considering Equation 1, even for comparatively low 
voltage levels, the current density passing through the 
insulator matrix is nonlinearly dependent on the electric 
field. For vanishing voltages, when eV~0, Equation 1 
can be simplified as [11] 



V exp 



'~h~ 



\JlmX I 



(2) 



The numerical evaluation of Equation 1 given in Figure 3 
shows the resistivity per unit area of polymer with 4-nm 
thickness and a quantum tunneling barrier height A of 
1 eV, with respect to the normalized voltage eV/A. This 
graph indicates that a polymer only exhibits close to 
ohmic behavior when subjected to low electric fields, that 
is, the resistivity of the polymer is approximately con- 
stant in a small region near the ordinate axis (see inset in 
Figure 3), permitting the use of the linear approximation 
provided by Equation 2. 

In this study, a rectangular potential barrier was as- 
sumed to model the electrical behavior of the tunneling 
resistor. Tunneling resistivity is numerically evaluated for 
A = 0.5 ev employing Equation 2 and illustrated in Figure 4. 
The tunneling resistance is drastically dependent on the 
insulator thickness, that is, tunneling resistance is sharply 
increasing as the insulator thickness is increasing. A cutoff 
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Figure 3 Polymer resistivity per unit area versus normalized 
voltage. The inset shows approximately ohmic behavior for low 
electric fields. 



distance can therefore be approximated at which tunnel- 
ing resistors with length greater than this threshold do not 
appreciably contribute toward the overall conductivity of 
the nanocomposite. In [12] and [13], the cutoff distance 
was assumed to be 1.0 and 1.4 nm, respectively. It is ex- 
pected that the resistivity of the insulator film is decreas- 
ing as the electrical field is increasing; so, when dealing 
with higher voltage levels, tunneling resistors with length 
greater than these cutoff distances may play a role in the 
nanocomposite conductivity. Hence, it was conservatively 
assumed in this study that tunneling resistors with length 
less than 4 nm contribute toward the nanocomposite 
conductivity. 

In the first step of this work, a three-dimensional con- 
tinuum percolation model based on Monte Carlo simu- 
lation was used to study the percolation behavior of 
an insulator matrix reinforced with conductive nanopla- 
telet fillers. Additional details on this modeling approach 



can be found in an earlier publication [14]. In the simu- 
lation, circular nanoplatelets are randomly generated 
and added to the RVE. The shortest distance between 
adjacent particles is calculated, and particles with dis- 
tance between them shorter than the cutoff distance 
are grouped into clusters. The formation of a cluster 
connecting two parallel faces of the RVE is considered 
the formation of a percolation network that allows 
electric current to pass through the RVE, rendering it 
conductive. 

Finite element modeling 

To study the electrical properties of nanocomposites, in 
particular their conductivity behavior, the employed 
modeling approach further involved the creation of a 
nonlinear three-dimensional finite element resistor net- 
work. Considering the excellent conductivity of the con- 
sidered nanoplatelets (e.g. tr= 10^ S/m for graphene), the 
electrical potential drop across the nanoplatelets was 
neglected. The tunneling resistors were modeled as non- 
ohmic, which resistivity is governed by Equation 1. 
Employing appropriate finite element formulations, the 
governing equation of an electrical resistor can be writ- 
ten as 



Vi 

v. 



(3) 



where /y is the electrical current passing between the 
ith and jth node; /Cy is the conductance of the resistor 
between nodes i and /; and Vi is the voltage of the ith 
node measured with respect to a node connected to 
ground. The system of the nonlinear equations govern- 
ing the electrical behavior of the nanocomposite was ob- 
tained by assembling the governing equations for the 
individual elements. The resulting nonlinear system of 
equations was solved employing an iterative method. 
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Figure 4 Tunneling resistivity versus insulator thickness. 



Results and discussion 

Modeling results 

The developed model was employed to investigate the elec- 
trical behavior of a polymer with A = 0.5 ev made conduct- 
ive through the uniform dispersion of conductive circular 
nanoplatelets with a diameter of 100 nm. In the simula- 
tions, the size of the RVE was chosen to be nine times the 
diameter of the nanodisks, which was ascertained to be 
large enough to minimize finite size effects. In an earlier 
study [15], the authors showed that the Monte Carlo simu- 
lation results are no longer appreciably RVE-size dependent 
when the RVE size is about eight times the sum of 2R + dt, 
where R and dt are the radius of the nanoplatelets and tun- 
neling distance, respectively. 

The graph in Figure 5 depicts the effect of filler loading 
on nanocomposite conductivity. As expected, a critical 
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volume fraction indicated by a sharp increase in nano- 
composite conductivity, i.e., the percolation threshold, can 
be inferred from the graph. 

In the following, electric current densities passing through 
the nanocomposite RVE were computed for different 
electric field levels and filler volume fractions. As il- 
lustrated by Figure 6, the current density versus volt- 
age curves were found to be nonlinear. The depicted 
electrical behavior of the conductive nanocomposite is thus 
clearly governed by the applied voltage in a nonohmic 



manner, which, as mentioned above, matches the expect- 
ation for a conductive nanocomposite at higher electric 
field levels. 

Figure 7 shows the variation of resistivity as a function 
of the applied electric field E in order to compare the 
nonohmic behavior for nanocomposites with different 
filler loadings. Note that resistivity values were normal- 
ized with respect to a reference resistivity measured at 
E = 0.8 V/cm. The results as displayed in Figure 7 indi- 
cate that the magnitude of the applied electric field plays 
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Figure 7 Normalized resistivity of nanocomposites with 100-nm nanodisks as a function of the applied electrical field. 



an important role in the conductivity of nanoplatelet- 
based nanocomposites. The employed modeling approach 
predicts nanocomposite resistivity to be a nonlinear func- 
tion of the applied voltage. Clearly, the nanocomposites 
exhibit nonohmic behavior where the resistivity de- 
creases with increasing voltage. Interestingly, the results 
in Figure 7 also indicate a reduced drop in resistivity and 
decreased nonohmic behavior for nanocomposites with 
higher filler volume fraction, that is, nanocomposites with 
higher filler loadings are less sensitive to the applied elec- 
tric field level. 



Comparison with experimental data 

To corroborate the simulation results, conductive epoxy 
nanocomposite samples were produced by in situ 
polymerization and their electrical behavior assessed 
as illustrated by Figure 8. Bisphenol-A epoxy resin and 
non-MDA polyamine curing agent (EPON 826 and 
EPIKURE 9551, by Hexion Specialty Chemicals, Columbus, 
Ohio, USA) were used for the fabrication of samples that 
were made electrically conductive by dispersing gra- 
phene nanoplatelets (xGnP-M-25, by XG Sciences, Lansing, 
Michigan, USA). 
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Figure 8 Normalized resistivity data versus applied electrical field from experiments with nanographene/epoxy samples. 
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Figure 9 Resistivity of nanocomposites with 100-nm circular 
nanoplatelets as a function of electric power. 



Graphene nanoplatelets were dispersed in acetone by 
sonication using a probe sonicater in an ice bath. In the 
following, epoxy was added to the mixture and sonic- 
ation was repeated. The solvent was evaporated by heat- 
ing the mixture on a magnetic stir plate and stirring 
with a Teflon-coated magnet. Remaining acetone was re- 
moved by using a vacuum chamber. The curing agent 
was added to the mixture and mixed with a high-speed 
mechanical shear mixer. The mixture was again degassed 



using the vacuum chamber and subsequentiy poured into 
a mold. A 2-h cure cycle was then performed at 120°C. 
Resulting samples were machined into circular disks with 
30-mm diameter and 3-mm thickness. The sample volume 
resistivities were measured at different applied voltages 
employing a Keithley 6517A electrometer connected to a 
Keithley test fixture (Keithley Instruments, Cleveland, 
Ohio, USA). 

Data in Figure 8 depicting the resistivity behavior of the 
epoxy nanocomposite samples was normalized with re- 
spect to the resistivity measured at an applied voltage of 
10 V. Samples with 1 and 1.25% graphene volume fraction 
exhibited high resistivity levels indicating a filler loading 
below the percolation threshold. For higher graphene vol- 
ume fractions of 1.75 and 2.25%, measurements indicated 
that percolation was achieved, and resistivity was found to 
decrease with the increase of the applied electric field. 
As predicted by the preceding modeling work, sample 
resistivity was found to be less sensitive to the ap- 
plied electrical field for higher filler loadings. Hence, 
modeling and simulation results are qualitatively in good 
agreement, indicating the validity of the assumptions 
undertaken for the numerical modeling. However, the data 
presented in Figures 7 and 8 also signify that further stud- 
ies are warranted to establish a quantitative agreement 
between numerical and experimental results. 

Characterization of resistivity behavior 

Gorrasi et al. [5] and Liu et al. [16] showed that the re- 
sistivity of carbon nanotube-based nanocomposites as a 
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function of the electric power P=Vx I can be described 
by an exponential expression: 

p = rP" (4) 

where a is an index which generally varies between - 1 
and 0. The value of a is indicative of the nonlinearity of 
the current-voltage relationship, i.e., a = 0 corresponds 
to ohmic behavior, and a decreases with increasing non- 
linearity of the current-voltage curve; r is a parameter 
relating to the resistivity of the nanocomposite when the 
electrical power passing through the sample is 1 W [16]. 

Computed nanocomposite resistivities are displayed as 
a function of the electric power in the graph in Figure 9. 
Data obeying Equation 4 appear in the form of straight 
lines owing to the graph's logarithmic scale. As shown in 
Figure 9, the slope of the lines decreases as the nonline- 
arity is decreasing with increasing filler loading. The 
values of a as a function of filler volume fraction are pro- 
vided in Figure 10. It is shown that a values are increasing 
with rising filler volume fraction. A discontinuity in a 
values can be observed in this graph for filler volume frac- 
tions of about 5%, which is associated with the percolation 
volume fraction. The behavior of data simulated herein is 
qualitatively congruent with results reported in [5] for car- 
bon nanotube nanocomposites. 

Conclusions 

In this study, the current-voltage behavior of conductive 
nanoplatelet-based nanocomposites was investigated. To 
this end, a numerical modeling approach was developed. 
The simulations predicted the resistivity of nanoplatelet- 
based nanocomposites to be strongly affected by the 
applied electric field. The nanocomposites exhibit non- 
ohmic behavior, that is, resistivity is a nonlinear function 
of the applied electric field. Further, nanocomposite re- 
sistivity was ascertained to decrease with increasing volt- 
age, while the degree of nonlinear behavior was found to 
decline with rising filler volume fraction. A good quali- 
tative agreement was observed between simulations 
and experimental data, the latter of which was obtained 
employing measurements on nanographene/epoxy nano- 
composites. The qualitative agreement between numerical 
and experimental studies encourages conducting a more 
comprehensive study to establish a quantitative agree- 
ment. The analysis further revealed that nanocomposite 
resistivity as a function of electrical power can be de- 
scribed by an exponential relation, where the exponent is 
a measure of the deviation from nonohmic behavior of the 
conductive nanocomposite. 

Nomenclature 

dt, tunneling distance 
e, electron charge 



E, applied electric field 

h, Planck's constant 

/, electrical current 

/, quantum tunneling current density 

k, conductance of tunneling resistor 

L, dimension of representative volume element 

m, electron mass 

P, electric power 

r, resistivity parameter 

R, radius of circular nanoplatelets 

Si, S2, limits of barrier at Fermi level 

As = S2 - 5i 

V, voltage across insulator 

filler volume fraction 
a, nonlinearity factor 
y6, correction factor 

A, height of barrier for quantum tunneling 
A , mean barrier height for quantum tunneling 
p, electrical resistivity 
a, electrical conductivity 
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